Aging in the linear harmonic oscillator 
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The low temperature Monte Carlo dynamics of an ensemble of linear harmonic oscillators shows 
some entropic barriers related to the difficulty of finding the directions in configurational space which 
decrease the energy. This mechanism is enough to observe some typical non-equilibrium features of 
glassy systems like activated-type behavior and aging in the correlation function and in the response 
function. Due to the absence of interactions the model only displays a one-step relaxation process. 
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Slow relaxation processes are widespread in condensed 
matter physics. These include magnetic relaxation in 
spin glasses, transport processes in structural glasses, 
pinning effects in superconductors among others. A large 
class of these systems show what is commonly referred as 
aging, i.e. dependence of the response of the system on 
the time in which it is perturbed. Aging effects |^] are a 
signature that the system is far from thermal equilibrium 
and consequently the fluctuation-dissipation theorem is 
not valid ||]. It has been realized quite recently that ag- 
ing is indeed a solution of the off-equilibrium dynamics in 
some exactly solvable models . Aging appears if re- 
laxation to the equilibrium is slow due to the presence of 
energy barriers in a rugged free energy landscape as well 
as in systems with entropy barriers with a quite simple 
landscape [||J|] . In this last the system relaxes to- 

wards the equilibrium, the number of directions in phase 
space where the system can move decreases progressively. 
This means that the system needs more time to decorre- 
late or to forget the previous configuration. This effect 
is usually encoded in the two time correlation function 
where the C{t,t') depends on both time indices |^]. 

From previous considerations it is clear that aging 
can also be present in extremely simple relaxing sys- 
tems without any interaction, the only condition being 
the progressive reduction of available phase space where 
the energy decreases. This was an essential ingredient 
in the Backgammon model recently proposed to explain 
glassy behavior in the absence of energy barriers |^ . Here 
we consider a simpler example and analyze the Brownian 
oscillator. The Brownian oscillator is usually studied in 
the Langevin approach. It is described in any textbook 
of stochastic theory It is possible to show that in 
this case there are no slow processes involved. In fact, 
the relaxation turns out to be exponential as expected for 
the dynamics of a particle in a single parabolic potential 
well. Here we consider the Monte Carlo approach and 
choose a dynamics based on the Metropolis algorithm M . 



This Monte Carlo approach was already studied in a dis- 
ordered model with long-range interactions which turns 
out to be non trivial, at least in the zero temperature 
limit [0 . The simplest case of an harmonic oscillator is 
solvable and we analyze the dynamics here. 

In we checked that, after a suitable rescaling of 
time, the equilibrium Langevin and Monte Carlo dynam- 
ics are equivalent. Also, we showed how the Langevin 
dynamics can be derived from the Monte Carlo dynam- 
ics in the limit of small changes. Here, we will see that 
the same results are valid. We will obtain the dynami- 
cal equations for the energy, the correlation and the re- 
sponse function. We will also study the low temperature 
dynamics, showing the similarities and differences with 
more realistic models for glasses. 

The harmonic oscillator has an energy. 



E = -Kx^ 
2 
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where K is the Hooke constant and x defines the posi- 
tion of the harmonic particle. Let us consider an ensem- 
ble of N independent linear oscillators with total energy 
E{{xi}) — {K/2)J2iXi- A change of {xi} is proposed 
{xi —> x[ = Xi-V ri/y/N,Wi} where {ri} is randomly cho- 
sen for each oscillator from a Gaussian distribution of 
zero mean and finite variance equal to A^. The change 
is accepted with probability 1 if the energy decreases, i.e. 
if SE = E({x'j^}) — E{{xi}) is negative. Otherwise the 
change is accepted with probability exp(— /3(5£') where 
/3 = ^ is the inverse of the temperature of the heat bath. 

Let us sketch the main derivation of the dynamical 
quantities jl^. We first consider the probability that a 
given set of movements {xi —>■ x'^ — Xi +ri/^/N} changes 
the energy in a quantity SE. This probability is given 
by. 
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For simplicity we have considered the case in which the 



mean position M = ■^5]^]^(a:i) of the initial condition is 
zero. The average (...) is done over different dynamical 
histories starting with the same initial condition for the 
ensemble. Using the integral representation for the delta 
function in the thermodynamic limit ^ oo we obtain, 



P{SE) = i'iTTKEA'^y^ exp( 
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Because the probability distribution P{SE) only de- 
pends on the energy itself the dynamics is then Marko- 
vian and simple to solve. Obviously this result is solely 
due to the simplicity of the model. According to the 
Metropolis dynamics the equation of evolution for the 
energy is, 



dE 
'dt ' 

which yields 



dxxP{x)+ / dx xP{x) exp{—Px) (4) 
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where a = {KA^/16E)2^ = ^ and 

evic{x) ^ / exp(-a;2)dx (6) 

J X 

fit) = ac^e-'^''^(i~2£(t)/3) e^ic{a{t){AE{t)l3 - 1)) (7) 

It is easy to check that the only stationary solution of 
this dynamical equation corresponds to the equilibrium 
solution with E = T/2 (in agreement with the equipar- 
tition theorem). The equation (||) is already closed and 
yields the dynamical evolution of the energy at all times. 

Knowing the evolution of the energy we can calculate 
the acceptation rate. This is defined by 
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In equilibrium, we have E^q — ^, and the acceptation 
rate becomes A^^g — erfc(aeg). 

Similarly one can derive equations for the correlation 
C(t,t') and response function G(<, i') (from now on, we 
will consider the first time t' as the smallest one, i.e. 
t' < t) defined by. 



Cit,t')^-{Y,x.it')x,it)) 
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where M(t) — ii/N)J2i{^i{t)) is the average position 
of the ensemble of oscillators. The response function is 
computed with the energy E = iK/2)^^xj — h^-Xi 
starting from an initial condition ilf = at zero field 
and taking finally the limit h 0. In the case that 
Mit — 0) ^ the computation is more involved and 
the equation of motion for the response function Git,t') 
involves also the energy response function Gsitjt') = 

( Sh{t'] ) ■ ^'^^ difficult to obtain the equations for 
both correlation and response functions. Proceeding in a 
similar way as for the case of the energy we get. 



dCit.t') 
dt 
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where ait) and fit) were previously defined. The dif- 
ference in the equations for the G and G concern only 
the initial condition. Note that again the Markovian 
properties of the dynamics are manifest because the time 
derivative of the G (or G) depends solely on the C(i,t') 
(or G(t, i')) itself and the initial conditions C(t',i') = 
2Eit')/K, limt^(^t')-Git,t')^0. 

We can easily integrate the equations for the correla- 
tion and response functions. We obtain that they depend 
on the energy (through the function fit)) at all previous 
times, 



Cit,t') = Slexp(- 
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With these exact results we can also calculate the 
fluctuation-dissipation parameter. 
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where a' ^ ait') and E' = Eit'). Note that X(i, t') only 
depends on the smallest time t' . 

It is straightforward to check that in thermal equilib- 
rium both correlation and response only depend on the 
difference of times, the fluctuation-dissipation theorem 
TGit) = -Git) is satisfied and the Xit,t') = 1. This is 
a general consequence of the detailed balance property 
inherent to the Metropolis algorithm. With an appro- 
priate rescaling of time we find also that the equilib- 
rium form of correlation and response function are the 
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same as in the Langevin case. Equations (|5|) and ( p_l[|l2| ) 
are much different that the corresponding ones in the 
Langevin dynamics. In particular, the response function 
at equal times is 1 in the Langevin dynamics while in the 
Monte Carlo case it is different from 1, even in thermal 
equilibrium. 

The Langevin case is a limit of the Monte Carlo dy- 
namics. This result has been demonstrated in the frame- 
work of the SK spherical spin-glass model and also 
applies here. When the size of the typical movement A 
goes to zero equations (p|jll|,p^) become, 
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with initial conditions C{t,t') = KE{t')/2 and 
limf^(^fi-^-G{t' ,t') = 0. These are the same equations as 
in the Langevin dynamics with uncorrelated white noise 
with variance 2T if the time is rescaled by the quantity 
A^/3/2. This means that both dynamics are essentially 
equivalent in case the rescaling factor is finite. 

The interesting case corresponds to the low temperature 
limit P oo for A fixed. In this case a new relaxational 
dynamics driven by a low acceptance rate is found. Lin- 
earization of eq. (||) around the equilibrium solution yields 
a relaxation time which diverges at low temperatures like 
r ~ /32 ex-p{KA'^f3/8). This implies a divergence of the 
relaxation time of an activated type. 
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FIG. 1. Cooling experiment. Values of the energy when 
we decrease and increase the temperature of the system at 
different ratios. From top to bottom, ratios 0.1,0.01,0.001 
and 0.0001. 



At very low temperatures the harmonic oscillator 
wants to relax to a configuration of very small entropy 
(indeed, because the oscillator is classical, the entropy 
diverges like log(T) at low temperatures). In this situa- 
tion the oscillator spends the major part of time looking 
for the ground state configuration. In some sense, the 
dynamics generates itself entropy barriers in a single po- 
tential well. This means that if we perform a cooling 
experiment, decreasing and increasing the temperature 
at a fixed rate, we expect that the system fails to relax 
to the equilibrium energy (see figure ^ . This is a typical 
feature of glassy systems. 

Let us consider the evolution of the energy at zero tem- 
perature. In this case, only those changes {Sxi} which 
decrease the energy are accepted. Close to the equilib- 
rium point X — the system will relax very slowly, mainly 
because the largest part of the movements are rejected. 
The relaxation of the energy at zero temperature is given 
by, 
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To obtain the long time behavior we expand the error 
function in the limit a — > oo. 
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We should note that previous equation is extremely 
similar to that derived in the Monte Carlo dynamics of 
the Sherrington-Kirkpatrick spherical model in the adi- 
abatic approximation at zero temperature. In terms of 
the parameter a (defined after eq.(|^)) the equation ( [20| ) 
can be written in the simple form. 



da exp(— a^) 
'dt " 



(21) 



For large times, the parameter a grows logarithmically 
in time. 



It ^ 2f 1 

a(t)^(log(— )-f-log(log(— ))^ 

the acceptation rate eq.(^) decays like 
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plus subdominant logarithmic corrections. The energy 
also decays logarithmically in time 
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Similarly the correlation function satisfies the equation 
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The same equation is fulfilled for the response func- 
tion. Using the asymptotic differential equations for the 
energy (pO|), the correlation ( p5| ) and the response func- 
tion, we can show that the correlation function for long 
times displays a solution of the type C{t, t') ~ ^tI^^j 



and G(t,t') ~ -^=]^g{t), being g{t) = exp(-a^). Using 
the asymptotic expression for the energy (|^ we get, 

' 2t' 
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FIG. 2. Correlation function for different waiting times. 
From top to bottom, = 3 ■ 10^,1 ■ 10^,1 • 10**,1 • 10^. The 
short lines show the calculated asymptotic behavior eq(p6|) 

This approximation is valid in the asymptotic limit 
of large values of t' . The normalized correlation function 
shows aging behavior with a simple scaling form t' /t plus 
some logarithmic corrections. Apparently the response 
function ( ^ ) does not show aging because it does not 
depend on a ratio of functions depending on t and t' . But 
this is an artifact of the normalization factor necessary 
to make the response function to take a finite value at 
equal times. In fact, the leading behavior of the response 
function decays to zero for large values of t and an ap- 
propriate normalization of the response function at equal 
times is necessary (in the same way as has been done for 
the correlation function). Note that for Langevin dy- 
namics the normalization of the response function is not 
necessary since the G{t, t') at equal times already takes 
a finite value by definition (e.g., limt/^(()- G{t,t') = 1). 
The normalized response function takes the simpler form, 
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G{t,t') 
G{t' + 0, t' 



^( 



(28) 



which displays aging with the same leading behavior as 
the normalized correlation function. We can obtain in- 
formation about the dynamics (and in particular, about 
the response function) from the remanent magnetization 
. In the present model the magnetization corresponds 
to the average position of the ensemble of oscillators (de - 
fined after (p^). The main equations (||),(|ll|) and ( |T^ ) 
have been derived in the absence of external field and 
starting with zero initial magnetization. In this case it is 
natural to compute the zero-field cooled magnetization. 
In this procedure the system is suddenly cooled down to 
a given temperature and after a waiting time a small 
step field h{t') = hQ{t' — i^) is applied. If the value 
of h is small enough then we are in the linear response 
regime. The magnetization starts to grow according to 
the relation, 



ZFC 



G{t,t')h{t')dt' = 



= h / G{t,t')dt' = hl{tyj,t) 



(29) 



The quantity I{tw,t) = G{t,t')dt' defines the in- 
tegrated response function. From the exact expression 
obtained for the response function (|lj) we get, for the 
integrated response function, 

I{tn.,t)^^{l-cxp{- f{t')dt')) (30) 

In the large time limit t ^ oo the zero-field cooled mag- 
netization converges to its equilibrium value, the field- 
cooled magnetization Mpc- It is easy to check, from 
d^), (|0|) that Mpc is given by Mpc = ie, the equi- 
librium linear magnetic susceptibility = Mpc/h = 
is independent of the temperature. 

Another interesting quantity to be calculated is the 
anomaly in the response function defined as. 



X= lim / dt'G{t,t') 

t^oo 



Geq{T)dT 
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Geqir) = lim G{tw,tn,+T) 



(31) 
(32) 

(33) 
(34) 



For a finite value of /3, the system decays to equilibrium 
in a finite time and for long times the integral dt' f{t') 
behaves as feqt. This implies that the magnetization re- 
laxes exponentially to zero, showing no aging for large 
values of and the 'anomaly' relaxes exponentially to 
zero too. The behavior of the anomaly and the zero-field 
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cooled magnetization is more interesting at zero tempera- 
ture. In that case it can be shown that the leading behav- 
ior of the anomaly decays algebraically (as l/Kt) to zero. 
Using the asymptotic behavior of the energy eq.(p^, it 
is easy to check that the zero-field cooled magnetization 
goes like, 

MzFC t,„/ log(7f) a 
Mfc ^ t Vlog(^)J 

Using the linear response relation Mzfc + Mtrm ~ 
Mpc where Mtrm is the thermo-remanent magnetiza- 
tion obtained by quenching the system in an (small) ap- 
plied field and removing it at t^^ we get 



(35) 



TR.M 



Mfc 



^ ^log(^) 



(36) 



Both Mzfc ^md Mtbm show aging with the leading 
t/tw scaling behavior. In figure 3 we show the thermo- 
remanent magnetization for the oscillator model for dif- 
ferent values of L„. 



FIG. 3. Thermo-remanent magnetization for different wait- 

From left to right, 



ing times obtained from eq.(R 
U = 10,10^10^10'*,10^ 



It has been suggested that the X {t,t') c ould be inter- 
preted as an effective temperature |L4| , [L5| . If we define 
Tf{t') = X{t,t')/T then, from cq. (|l|), the fluctuation- 
dissipation theorem is obeyed with the effective temper- 
ature Tf{t'). While this is a formal relation it would 
be interesting if the effective temperature derived in this 
way had some deep physical meaning. On the other 
hand, a well-founded physical interpretation of the viola- 
tion of the fluctuation-dissipation relation, to our knowl- 
edge, does not exist. Note that it is possible to de- 
fine different fluctuation-dissipation ratios (for instance. 



X{t,t') = -TG/^) all giving X{t,t') = 1 in equilib- 
rium. The definition here adopted is the conventional 
one which allows to obtain a closed expression for the in- 
tegrated response function in case the function X is solely 
function of the correlation C(t, t') 1^,1). On general phys- 
ical grounds one would expect an effective temperature 
larger than the temperature of the bath. To raise the 
temperature should contribute (by the equipartition the- 
orem) those degrees of freedom which, during the process 
of relaxation towards the equilibrium, still are not frozen. 
For the simple model considered here such an interpreta- 
tion seems to work. From equation (^5|) it can be shown 
that the effective temperature for a system relaxing at 
zero temperature is given by the relation Tf{t') 2E{t'). 
Consequently the effective temperature and the dynam- 
ical energy in the off-equilibrium regime are related by 
the thermodynamic relation suggesting that some kind of 
adiabatic theorem holds for this simple system in the long 
time limit. One can then ask if the whole time dependent 
probability distribution pt{{xi}) in the long-time limit is 
of the Boltzmann type but dependent on an effective tem- 
perature Tf{t), i.e. pt{{x^}) ~ exp(-£:({x,})/r/(t)). It 
is easy to check that such a result is not possible [0 and 
equipartitioning is valid only for some finite moments of 
the probability distribution (for instance the second mo- 
ment, i.e. the energy). 

In conclusion, we have studied the Monte Carlo dy- 
namics of an ensemble of linear harmonic oscillators. The 
extreme simplicity of this model makes it exactly solv- 
able without loosing the interesting features of the non- 
equilibrium dynamics driven by entropic barriers. In this 
way, we are able to gather quite a lot of information and 
derive all relevant dynamical quantities with reasonable 
analytical effort. 

We find a very slow relaxation near zero temperature, 
driven by a low acceptance rate, similar to that found in 
the Backgammon model [|j , models of adsorption jl| 1 a s 
well as models for compaction of dry granular media 1 19| . 
In these cases, the origin of the slow relaxation is the ex- 
istence of entropic barriers, although they are set up by 
different mechanisms. Note that the notion of entropic 
barrier or entropic trap is quite similar to the concept of 
effective volume in free volume theories. In our case this 
manifests as a inverse logarithmic law decay of the en- 
ergy eq.(p4|) while in compaction of granular media this 
decay is found for the density of compaction. The model 
has also in common with models for glasses aging in the 
correlation function for long times. The correlation func- 
tion C{t, t') presents a j behavior with some logarithmic 
corrections (with t' the smallest time in the correlation 
function). It is interesting to note that this corrections 
appear also in the Backgammon model (and presum- 
ably also in adsorption models |l8| and models for com- 
paction of dry granular media |1£ ) but do not appear in 
models with Langevin dynamics |20 21 1). We have found 
also aging in the magnetization (the integrated response 
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function). This behavior appears associated to the alge- 
braical decay of the 'anomaly' as t goes to infinity. Due 
to the zero value of the anomaly we expect a finite value 
of the overlap between two replicas in the large t 
limit if putted in the same configuration at (cloning 
procedure). This expectation stems from the simplicity 
of the landscape in this model (a single parabolic well). 
Consequently, this model falls into the first dynamical 
category (class I) proposed in Q. However this model 
shares some features of the p-spin model (with p > 2, be- 
longing to class II) like aging in the integrated response 
function. Furthermore, this simple model lacks a fast 
process decaying to a plateau and also a two time de- 
pendence of the fluctuation-dissipation parameter, which 
only depends on the smallest time. This is probably due 
to the absence of a cooperative behavior. 

Acknowledgments. We are grateful to S. Franz 
and Th. M. Nieuwenhuizen for a careful reading of the 
manuscript. The work by L.L.B. and F.G.P. has been 
supported by the DGES of Spain under grant PB95-0296. 
The work by F.R has been supported by FOM under con- 
tract FOM-67596 (The Netherlands). 



ence on glassy systems, June 1996, edited by M. Rubf, 

Preprint cond-mat 9607224; 
[12] S. Franz and M. Mezard, Physica A 210 48 1994; 
[13] A. Barrat, R. Burioni and M. Mezard, J. Phys. A (Math. 

Gen.) 29, 1311 (1996); A. Barrat, Ph. D. Thesis (Uni- 

versite Paris VI), 1996. 
[14] S. Franz and F. Ritort, J. Phys. A (Math. Gen.) 30, L357 

(1997). 

[15] L. F. Cugliandolo, J. Kurchan and L. Peliti, Phys. Rev 

E55, 3898 (1997). 
[16] C. Godreche and J. M. Luck, J. Phys. A (Math. Gen.) 

29 1915 (1996) and Preprint cond-mat 9707052 
[17] This result can be obtained for Langevin dynamics, S. 

Franz and F. Ritort (unpublished) 
[18] J. W. Evans, Rev. Mod. Phys. 65, 1281 (1993) 
[19] E. Caglioti, V. Loreto, H. J. Herrmann and M. Nicodemi, 

A "Tetris"-like model for the Compaction of dry granular 

media Preprint cond-mat 9705195. 
[20] S. Ciuchi and F. De Pascuale, Nucl. Phys. B 300, 31 

(1988); L. F. Cugliandolo and D. S. Dean, J. Phys. A 

(Math. Gen.) 28, 213 (1995). 
[21] F.G. Padilla and F. Ritort, Langevin Dynamics of the 

Lebowitz-Percus Model, Preprint cond-mat 9703095 to 

appear in J. Phys. A. 



[1] L. Lundgren, P. Svedlindh, P. Nordblad and O. Beck- 
man, Phys. Rev. Lett. 51, 911 (1983); E. Vincent, J. 
Hammann and M. Ocio, in 'Recent progress in Random 
Magnets', ed. D. H. Ryan (Singapore: World Scientific) 
(1992), and references therein; L. C. E. Struik, Physical 
aging m amorphous polymers and other materials (Else- 
vier, Houston, 1978) 

[2] J. P. Bouchaud, L. F. Cugliandolo, J. Kurchan and M. 
Mezard, Out of Equilibrium dynamics in Spin Glasses 
and other Glassy systems Preprint cond-mat 9702070; 

[3] L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71 
173 (1993); J. Phys. A (Math. Gen.) 27 5649 (1994); 

[4] S. Franz and M. Mezard, Europhys. Lett. 26 209 (1994); 
Physica A 209 1 (1994); 

[5] F. Ritort, Phys. Rev. Lett 75, 1190 (1995); S. Franz and 
F. Ritort, Europhys. Lett. 31, 507 (1995); C. Godreche, 
J. P. Bouchaud and M. Mezard, J. Phys. A (Math. 
Gen.)28 L603 (1995) 

[6] A. Barrat and M. Mezard, J. Physique I (France) 5 

(1995) 941. 

[7] H. Rieger, J. Phys. A (Math. Gen.) L615 (1993); 

[8] N. G. Van Kampen, Stochastical Process in Physics and 

Ghemistry North Holland (1981); 
[9] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. 

H. Teller and E. Teller, J. Chem. Phys. 21 1087 (1953); 
[10] L.L. Bonilla, F.G. Padilla, G. Parisi and F. Ritort Phys. 

Rev. B 54,6 , 4170 (1996); Europhys. Lett. 34 3 159 

(1996) ; 

[11] E. Vincent, J. Hammann, M. Ocio, J. -P. Bouchaud and 
L. F. Cugliandolo, in Proceeding of the Sitges Confer- 



6 



